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Among transport-limited reactions, reactions involving polymeric chains play an important role. 
Both intramolecular reactions such as cyclization and intermolecular reactions have been extensively 
studied experimentally and theoretically, and have been shown to lead to complex kinetics. Despite 
these considerable efforts, there is to date no exact explicit analytical treatment of transport- 
limited polymer reaction kinetics, even in the case of the simplest model of flexible polymer - a 
phantom Rouse chain of monomers connected by linear springs. The main difficulty arises from 
the fact the motion of a single monomer in the chain is non Markovian. Here, we introduce a new 
analytical approach to calculate the mean reaction time of polymer reactions that encompasses the 
non Markovian dynamics of the problem. A key step of our method relies on the determination 
of the statistics of the polymer conformation at the very instant of reaction, which provides as a 
by product new information on the reaction path. We show that the typical reactive conformation 
of the polymer is more extended than the equilibrium conformation, which leads to reaction times 
significantly shorter than predicted by existing Markovian theories. Together, these results provide 
a better understanding of the complex kinetics of polymer reactions involved for example in the 
formation of loops of RNA or polypeptides chains. 

Reactions involving macromolecules and in particular polymer chains are ubiquitous. An important example of 
intramolecular polymer reaction is provided by cyclization reactions, which consist in forming a loop joining the two 
ends of a polymer. Such reactions have been widely studied both theoretically [IHH] and experimentally [51116] . mainly 
because of their relevance to biological processes. Indeed, the formation of loops and hairpins in DNA is a key process 
in the regulation of gene expression [17] ; in the context of protein folding, the cyclization of a polypeptide chain can 
be seen as an elementary step of the folding pathway |10) . Among intermolecular polymer reactions, search processes 
involving a polymer chain and a given target, be it a catalytic site or a pore in a confining cavity, play a prominent 
role, as exemplified by gene delivery or viral infection, which involve a step that is kinetically limited by the search 
for a nuclear pore by a nucleic acid [X8H50] . 

The theoretical description of polymer reaction kinetics requires to take into account the intrinsically complex 
dynamics of a polymer chain. The motion of a monomer depends on the dynamics of the entire chain, and therefore 
cannot be described as a Markov process. This non Markovian feature induces the emergence of multiple time scales, 
and can lead to subdiffusion [5T] [55] and non trivial reaction kinetics [7J [531 HI] ■ Numerous studies have been devoted 
to the theoretical analysis of polymer reaction kinetics, but until now all available explicit results rely on Markovian 
approximations of this non Markovian problem. 

The benchmark theory in the field has been developed by Wilemski and Fixman [H [2] , and assumes that all the 
hidden degrees of freedom of the polymer reach instantaneously their equilibrium distribution. Another classical 
theory is the harmonic spring model, where the whole polymer chain is modeled by a single spring, with an effective 
stiffness that takes into account the entropic stiffness of the chain [3j [25l [26] . A third theoretical approach is based 
on the renormalization group theory and leads for infinitely long chains to perturbative results in e = 4 — d, where 
d is the space dimension - More recent contributions include a formal iterative solution in dimension 1 [57], or 
a more refined treatment of the correlations 128:. These theories capture some features of the anomalous polymer 
dynamics, but explicit results invariably make use of a Markovian approximation. This assumption inevitably leads 
to a restricted range of applicability of these approaches, and clearly fails as soon as the reaction time is of the same 
order as the polymer relaxation time. In fact, results of numerical simulations have proved to significantly differ from 
the available theoretical predictions in a broad range of parameters [3 155H5T] . 

Here, we propose a new approach that directly deals with the non Markovian character of the problem. The key 
step of our method relies on the determination of the statistics of the polymer conformation at the very instant of 
reaction, which has been disregarded so far. We show that typical reactive conformations are in marked contrast 
with the equilibrium conformations as is implicitly assumed in existing Markovian theories. Our analytical approach 
provides a very accurate determination of the mean reaction time for both intra and inter molecular reactions valid 
for any range of parameters, which significantly improves standard polymer reaction kinetics theories. These results 
open new perspectives in the understanding of the complex kinetics of polymer reactions. 




FIG. 1: Sketch of the examples of intramolecular and intermolecular reactions studied in this paper. What is the 
influence of the presence of various monomers on the reaction kinetics between reactants attached to particular monomers ? 
We address this question in this paper by considering the two examples of the cyclization reaction (a) and the reaction between 
a reactant attached to the end of a polymer and a fixed target in a confining volume (b). 

Results 

We consider the classical model of a Rouse chain of N monomers connected by linear springs of stiffness k. The 
monomers experience a frictional drag of coefficient £ and diffuse with a diffusion coefficient D = kgT/£ in the force- 
field created by their neighbors. Even if this minimal model neglects both hydrodynamic interactions and excluded 
volume effects, it captures the main features of polymer dynamics [51] [55]. We denote the microscopic time scale by 
t = Q/k, which is the typical relaxation time of a bond in the polymer, and the microscopic length by Iq = y^kgT/k, 
which is the typical length of a bond. We introduce the positions r,, i € {1, N} of the N monomers, where quantities 
in bold stand for vectors in the d-dimensional space. The probability P(r±, ...,r^,t) to find the polymer chain in a 
given configuration at time t satisfies the Fokker-Planck equation [211 E2] : 

N 1 N 
i=l * i=l 

where V, = d/dri and Fj is the force acting on the i th monomer. This force is related to the monomer positions 
by Fj = fc(rj + i — 2ri + r^_i) with the convention that i"o = i"i and r« + i = r^. This defines the matrix M such 
that Fi — ~k'Y^ =1 MijYj. The matrix M is tridiagonal positive symmetric, its eigenvalues are therefore positive and 
given by Xj — 2[1 — cos((j — l)ir/N)] with j G {1, ...,N}. The first eigenvalue is Ai = and is associated with the 
motion of the polymer center of mass. The largest relaxation time of the internal conformations of the chain is named 
the Rouse time t r = C/(fcA 2 ) = N 2 (/(kn 2 ). 

To cover both cases of intra and intermolecular reactions, we focus on two examples of reactions sketched in Fig. 
[T] that are associated to two observables R bs- First, the position of the first monomer R bs = 1*1, relevant to 
intermolecular reactions, which diffuses at large times with the same diffusion coefficient Dqm — D/N as the polymer 
center-of-mass. Second, the end-to-end vector R bs = r JV — r ij relevant to intramolecular reactions, which is not 
diffusive since the mean squared displacement relaxes in finite time to an equilibrium value. To quantify the reaction 
time for intra and intermolecular reactions, we aim at calculating the mean first passage time [32j for the variable 
Robs to reach a value ||R bs|| = a given an initial probability distribution of the polymer position Pi n i(|r)), where 
a is the reaction radius (see Fig. [l]). Here |r) represents the vector with N components (r 1; ...,rjv) that defines a 
conformation of the polymer. We will consider the cases where the initial distribution is the stationary distribution 
-PiniQr)) = P s tat(|r)), or the stationary distribution restricted to values of |r) such that ||R bs|! = -R^bs- tne case °^ 
intermolecular reactions, we introduce a large confining volume V in which the reaction takes place (see Fig. [TJa). 

While the dynamics of R bs is non Markovian, the evolution of the full polymer conformation |r) is Markovian and 
obeys a renewal equation which is the starting point of our analysis. Let us define /(|r),t) the probability density 
that, starting from the initial distribution, the reactive region is reached for the first time at t with a configuration |r). 
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The renewal equation then takes the following form which is valid for all the configurations |r) such that ||R b s || < a: 



P(|r),t|ini,0) = dt' d\r')f(\r'},t')P(\r},t-t'\ |r'),0). (2) 



Here, d\r) — dri...drjv, P(|r),t|ini, 0) is the probability of a configuration |r) at t in the absence of target when 
the initial distribution at t = is Pi n ;(|r)), and P(\r),t — t'\ |r'),0) is the probability of observing the configuration 
|r) at t given that the configuration |r') was observed at t = 0. We introduce the splitting probability distribution 
7r(|r)) = 7r(i*i, rw) that represents the probability density of observing a configuration |r) when the reaction takes 
place. Taking the Laplace transform of the renewal equation ([2| and developing for small values of the Laplace 
variable yields an integral equation that links the mean first- passage time r and the splitting probability 7r(|r)): 

/>oo 

rP stat (|r)) = / d*[P(|r>,i|7r,0)-P(|r>,i|ini,0)]. (3) 
Jo 

Here, we have introduced P(|r), t\ir, 0) the probability of a configuration |r) at t given that the configuration at t = 
is taken from the splitting probability 7r, which reads: 

P(\r),t\n,0) = J d\r')n(\r'))P(\r),t\ |r'),0). (4) 



The equations (3]4) together with the normalization condition for 7r(|r)) form an integral equation that completely 



defines ir and r, but which is very difficult to solve in the general case. Let us introduce a final position for the 
observable R obs that is located inside the reactive zone (||R obs || < a). Integrating Eq. ^ over the |r) such that 
Robs = R b s gi yes an exact expression of the mean reaction time r as a function of the splitting probabilities: 



rP stat (R obs )= / dt P(R f obs ,t|7r,0)-P(R obs ,i|ini,Q) . (5) 
Jo L J 



This exact expression generalizes the results obtained in Refs. [33H35] for Markovian systems. The term P s tat(R obs ) 
represents the stationary probability distribution of observing the observable with the value R obs - In the case of 
intramolecular reactions, this term is sometimes called the j- factor jT7j. In the case of intermolecular reactions, for 
which R obs = P sta t(R obs ) = l/V', where V is the confinement volume, and we will derive below the large volume 
asymptotics of the mean first-passage time. Note that, by construction, the formula ^ provides the same value of r 
for all possible values of R obs inside the reactive zone. 

The equation ^ shows that the calculation of the mean first-passage time requires the determination of the dis- 
tribution of the polymer configuration at reaction n, which is highly non trivial. A local equilibrium assumption, 
which turns out to give the same results as the Wilemski-Fixman approximation p] [21 129) , then consists in approx- 
imating the splitting probability by the stationary probability restricted to configurations |r) such that R obs lies on 
the surface of the target [7r(|r)) ~ P s tat(|r)| ||R bs|| = «)]• This approximation is Markovian because it assumes that 
all the variables relax instantaneously to an equilibrium distribution. Here we go beyond this Markovian assumption 
and keep track of non Markovian aspects of the problem by calculating n in a self-consistent way. As we proceed to 
show, this distribution tt of the polymer conformation at the instant of reaction markedly differs from the equilibrium 
distribution, showing that the non Markovian features of the kinetics cannot be ignored. For clarity, we present the 
method in dimension 1, and discuss the results in dimensions 1 and 3. 

In the 1-dimensional case one can take the target size a = and cc obs = without loss of generality. The 
key assumption of our approach is that the splitting probability tt(\x}) — tt(xi, x^) is a multivariate gaussian 
distribution that is fully characterized by the averages (which we call mf) and the covariance matrix (denoted a^A 
of the variables Xi. The moments of 7r can be calculated by using self-consistent equations, that are derived from ([3]) 
in the supplementary information (SI): 



oo 

P(0,t|7r,0)^'°- J P(0,t|ini,0)/i i ; ni, ° =0 (6) 

n/n j-l n\ ( ^ '•* i ^-0 stat.*\ r>/n xl* ■ n\ ( ini,* . ini.O ini.O stat,*\ n /«\ 

P(0,t|7r,0) (7^.' +Mi J - nM ini >°) (7ii +Mi ' Hj -o-y JJ =0 ( 7 ) 



dt 

OO 

dt 



o 



Here, /it^' and 7^'* are the moments of the distribution P(\x), i|a; obs = 0, t; tt, 0) (the distribution of the configurations 
\x) at t such that x obs = and given that the initial configuration is taken from the splitting probability). These 
moments /x"' and 7^'* can be related to the initial moments mf and afj by projection and propagation formulas 
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that are given in SI. In the same way, and 7™'* are the moments of the distribution P(\x),t\x b s — 0,f;ini, 0), 

whereas <7^* at '* is the covariance matrix of the distribution P s tat(|^)|^obs = 0). Last, the propagator P(0, t\ir, 0) can be 
written as a function of and crj (see SI). Equations ( 6|7 | provide a closed system of explicit equations that fully 
determines the unknown moments and tr^-, and therefore the distribution tt(\x)). These equations can be extended 
to the case of a 3-dimensional space (see SI). Together with Eq. ([5]), Eqs.(6|7) finally enable the determination of the 
mean reaction time and constitute one of the main results of our work. 



Discussion 

In order to test the validity of our non Markovian approach, we compared its predictions with the results of several 
numerical simulations in 1 and 3 dimensions for both cases of inter and intramolecular reactions (see Figs 2|3|4 1. 
The Brownian dynamics simulations were carried out by using the algorithms presented in Refs. |291 136) (see SI for 
details). Since the Wilemski-Fixman approach is the Markovian theory that matches best simulation results until 
now [29] . we also compare our predictions to this theory, which will be referred to below as the Markovian theory. 



O 




FIG. 2: Theory and simulations of the cyclization reaction. Comparison between the markovian theory (dashed lines), 
the non-markovian theory (continuous lines) and simulation results (symbols). Different colors correspond to different sizes of 
reactive region (blue: a/lo = v3/10, red: a/lo = v3, green: a/lo = 5y3). The squares and diamonds symbols correspond to 
the simulations data of Ref. [29] , whereas the circles correspond to our Brownian dynamics simulations. The non-markovian 
results are calculated by approximating the covariance matrix of the splitting probability by its stationary value. The units of 
length and time are Iq and £/fc. 



In the case of intramolecular reactions exemplified by the cyclization reaction (see Fig. [jj,), our non Markovian 
theory is in excellent agreement with the simulations for all sizes a of the reactive region and for all polymer lengths, 
and significantly improves the results of the Markovian theory (Fig[2 ). To our knowledge, our non Markovian approach 
is the first theory that provides an accurate description of the eye ization time for all ranges of parameters. In the 
case of intermolecular reactions where for example the reactive site is the first monomer of a chain (see Fig. ^p), our 
theory yields an excellent quantitative determination of the reaction time and significantly improves the results of 
existing Markovian theories (see Figs. [3^,j4^). As stated above, the key element of our analysis is the determination 
of the distribution of the polymer conformation at the very instant of reaction. The very precise determination of 
the mean first passage time is a consequence of the fact that the Gaussian approximation accurately describes the 
distribution n (see Figs[3]3,c and Figs|4jD,c). We stress that the Gaussian prediction of our non Markovian approach 
markedly differs from the equilibrium distribution implicitly assumed in Markovian theories. Remarkably, the reactive 
conformation of the polymer for both inter and intramolecular reactions is actually significantly more extended than 
the equilibrium conformation, which yields reaction times notably shorter than predicted by Markovian theories. 

We now derive analytical formulas for the mean reaction time from Eqs. ( 5|6|7 | in different limiting regimes in the 
most relevant case of dimension 3. These simple explicit formulas enable the definition of the different regimes of 
polymer reaction kinetics, and establish clearly the validity domains of the Markovian theories. Scaling relations rely 
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FIG. 3: Theory and simulations of a diffusive variable in one dimension: reaction of the first monomer of a 
Rouse chain with N = 10 monomers with a target located at x = in a confining volume V . (a) Mean first passage 
time (rescaled by the volume) as a function of the initial distance x\. Red symbols: results of numerical simulations for a volume 
V = 316io, error-bars are 95% confidence intervals. Lines: theoretical estimations of t/V (green: markovian approximation 
; red: non markovian theory), (b) Marginal splitting probability histogram for the 6 th monomer position (ir(xfi)) when the 
starting position is x\ — 5.7Zq . Lines: theoretical distribution predicted by the non Markovian theory (red) and the Markovian 
approximation (green), (c) Average position of the monomers at the reaction for x\ = 5.7Zo (simulation and theory, with the 
same colors and symbols as in (a)). The units of length and time are io and (/k. 



in part on the behavior of the mean squared displacement of a monomer with time: 

{6Dt iH<CA 
at 1 ' 2 if (/k < t < t r ~ N 2 C/k (8) 
6D CM t if t > t r 

where a is a numerical coefficient. The anomalous diffusion at intermediate time scales involves all the time scales of 
the Rouse chain and becomes important in the limit of long chains N ^> 1. At longer time scales, the behavior of the 
monomer is diffusive with diffusion coefficient Dcm- 

Small target regime a <C h/ \fN. In this limit of a target much smaller than the bond length, the reaction kinetics 
depends essentially on the short time properties of the search, and the mean reaction time can be shown from Eqs. 
( 5|6|7 l to asymptotically follow: 



{V/ (4irDa) (intermolecular reaction) 

•\/7r/8(ZoiV) 3 / 2 / (Da) (intramolecular reaction) . 

This shows in particular that for large N the mean cyclization time scales asr~ N 3 / 2 /a for small targets. The strong 
dependance of the mean reaction time with the target size in this limit is the signature of a non-compact exploration 
[7]. It is notable that in this non-compact limiting case, the non Markovian theory predicts the same result as the 
two classical Markovian theories (the Wilemski-Fixman theory [HE] and the harmonic spring model [25j [29] ) in 
the limit of small reaction radius, which validates Markovian approaches in this regime. 

Intermediate target regime lo/s/N <C a < loVN. In the regime of long polymer chains, a single monomer displays 
at intermediate time scales a subdiffusive behavior (Ar 2 ) ~ t 2 / dw [Eq. (Isb] , thereby defining a walk dimension d w = 4 
that is larger than the spatial dimension d — 3. As a consequence, a monomer is able to densely explore the space, 
and the time to reach a target much smaller than the polymer size Iq^/N is asymptotically independent of the target 
size a. In the case of intermolecular reactions, we show from Eqs. ( 5|6|7 ) that the reaction time averaged over all 
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FIG. 4: Theory and simulations of a diffusive variable in 3 dimension: reaction of the first monomer of a 
Rouse chain with N = 10 monomers with a target of size a = 1.7Zo centered around the position r = in a 
confining spherical volume V of radius 28.5Zo- The color code is the same as in Fig [3] (a) Mean first reaction time r(r?) 
(rescaled by the confining volume) as a function of the initial distance between the reactants. Symbols: simulations ; Lines: 
theoretical estimations (green: markovian approximation ; red: non-markovian theory). Dashed lines give estimations of r by 
using 1 1 Robs 1 1 = a i whereas the continuous lines are calculated with R obs = 0. All calculations assume the spatial isotropy of 
the covariance matrix, (b) Bars: histogram of the radial component of the position of the 6 th monomer at the reaction when 
r? = 28Zo- Lines: theoretical distribution predicted by the non Markovian theory (red) and the Markovian approximation 
(green), (c) Radial average positions of the monomers at the reaction when r° = 28Zo (red symbols), compared with the 
theoretical predictions of the markovian theory (green line) and of the non-markovian theory (red line). The units of length 
and time are lo and Q/k. 



initial conditions reads: 

V 



A_irD CM a c s{a,N) 



(10) 



where a e s(a,N) is an effective target size, which is of the order of the polymer size (a c g ~ ZoviV) and does not 
depend on the real size of the target a for a — > 0. This equation has a clear interpretation : the polymer only needs 
to approach the target at a distance comparable to Iq\/N \ and then the reaction takes place instantaneously due to 
the compact search at small length scales. Importantly, if both the early analysis of De Gennes [7] and the Markovian 



theory predict the same functional form (10), they fail to predict a correct estimate of the effective target size a e g, 
which is found in realistic regimes to be underestimated by a factor of two as compared to the non Markovian result 
(see SI). This difference can be understood from the fact that the polymer is much more extended at the instant of 
reaction than in its equilibrium conformation (Fig. [5]a). In the non Markovian description, the polymer center-of- 
mass therefore needs to approach the target less closely than in the Markovian theory, leading to a faster reaction 
kinetics. Similarly, we find in the case of intramolecular reactions that the mean cyclization time (averaged over 
stationary initial configurations) is given by 



c 



a \ N 2 C , , 



where again the numerical function c(a/^ov N) is typically underestimated by a factor of two by the Markovian theory 
(see SI). 

The faster kinetics predicted by the non Markovian theory for both intra and intermolecular reactions is a direct 
consequence of the non trivial out of equilibrium distribution of the polymer conformation at the instant of reaction (see 
Fig. [5]). While the Markovian theory implicitly assumes an equilibrium conformation, the non Markovian theory shows 
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(d) Conformation at reaction 
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FIG. 5: Predicted polymer conformations at reaction. Left: Average radial position of the monomers when the reaction 
takes place for: (a) the reaction between the first monomer and a target and: (b) the cyclization reaction (Continuous line: 
prediction of non Markovian theory, dashed line: Markovian approximation). For both reactions, we also plot the sketch of the 
polymer shape when the reaction takes place, and an example of conformation drawn from the splitting probability distribution 
(left) which is in marked contrast with the stationary distribution (right). The reaction is assumed to take place along the 
vertical axis. The position of a monomer in the chain is represented by a color code. The sketch of the polymer shape for 
the cyclization reaction is artificially extended in the horizontal direction for clarity. Parameter values: (a): N — 800 and 
a = 0.32/0^; (b): N = 800 and a = 0.094l ^fN. The unit of length is l . 



that the reaction takes place when the polymer is in fact much more extended than in its equilibrium conformation, 
thus increasing the effective reaction radius, as seen in Fig. [5jb. This contribution of non Markovian effects turns out 
to be quantitatively important, since in this regime of intermediate targets the Markovian approximation leads to an 
error in the estimate of the reaction time of roughly 100%. 

While we have here focused on the Rouse model, we stress that the fact that such non Markovian effects are 
characterized by non equilibrium polymer conformations at the instant of the reaction holds true for more general 
models of polymer dynamics. Our method can in fact be extended to general gaussian models, which play a key role 
in polymer dynamics and enable the modeling of various physico-chemical conditions. Examples of such gaussian 
theories include the "pre-averaging" approach of hydrodynamic interactions |22j . the approximate Rouse modes in 
the case of self-avoiding polymers |37j . and the description of semi-flexible chains and branched polymers with a 
gaussian theory J3SJ 02] ■ Qualitatively, we expect the non Markovian effects to be significant when the search at 
small time scales is compact, so that the transport step plays a crucial role in the kinetics. This includes the case 
of self-avoiding chains as well as chains with hydrodynamic interactions in theta solvent. Last, at the experimental 
level, the formation of hairpins in nucleic acids [9l [11] or the folding of polypeptide chains [10) constitute important 
examples of cyclization reactions. So far, only Markovian theories of such reactions have been used to interpret 
observations |101 111) and we anticipate that taking into account non Markovian effects as quantified by our approach 
could improve the quantitative analysis of experimental data. 

To conclude, we proposed a new theory of polymer reaction kinetics that takes into account the non Markovian effects 
that control the dynamics of polymers. This non Markovian theory gives results that are in quantitative agreement 
with numerical simulations for all ranges of parameters, and therefore significantly outperforms existing Markovian 
approaches. Our analysis reveals that the non equilibrated conformation of the polymer at the instant of the reaction 
has an important impact on the reaction kinetics. We show quantitatively that the typical reactive conformation 
of the polymer is more extended than the equilibrium conformation, leading to reaction times that are significantly 
shorter than predicted by existing Markovian theories. Together, our results provide a better understanding of the 
complex kinetics of polymer reactions involved for example in the formation of loops of RNA or polypeptides chains. 



Supplementary information 



Appendix A: Derivation of the self-consistent equations (5,6,7) of the main text. 
1. Definition of the Rouse modes and choice of units 



Before describing how to obtain the self-consistent equations (5,6,7) of the main text, we introduce the notion of the 
Rouse modes, that considerably simplifies the Fokker-Planck equation (1). We remind that the matrix My that links 

the forces on the monomers Pi to the positions v.j is defined by the relation: £F S ; = — fc^jLi ^ij r j- M is therefore 
the following N x N tridiagonal matrix: 



/ 1 -1 



M 



-1 2 



-1 



\ 



-1 2 



-1 



V 



0-12-10 
.. 0-12 -1 
.. .. 0-11/ 



(Al) 



Because M is symmetric positive, it can be diagonalized. The eigenvalues of M are: 

A,- = 2{1 - cos[(i - l)ic/N]} (A2) 

We can write M = QDQ -1 , with D the diagonal matrix with (Ai, Ajv) on the diagonal and the passage matrix Q 
is: 

(A3) 



Qn = 1/VN ; Q y = d -cos \{i - 1/2) (j - 1)tt/N] if j > 2 



N 

Note that the inverse of Q is its transpose matrix. We define the Rouse modes (ai, ...,ajy) by the relations: 

N N 



(A4) 



Note that the transformation ( |A4[ ) is equivalent to the Fourier transform when N is large. We study the observables 
R-obs that can be expressed as a linear combination of the positions r i; or equivalently of the Rouse modes a.j. An 
observable R bs is therefore associated with a set of N coefficients (bi, 6jv): 



■obs 



R 



N 



hen = (6|a) 



(A5) 



»=i 



In order to have more contracted notations, we omit the subscript of R bs and simply call the observable R (we reserve 
capital letter for the observable). In equation (A5), we have introduced the notation that |tt) is a column vector with 
iV components (ui, ...,Un), (u\ is its transpose, and (u\v) = X)i=i u i v i is the scalar product between the vectors \u) 
and \v) . Note that quantities in bold represent vector in the physical 3-dimensional space, to be distinguished from 
the N components vectors noted \u). If the observable is the end-to-end vector (R = rjv — ri), then the coefficients 
are hi = Qni — Qu- In the case that the observable is the position of the first monomer (R = ri), the coefficients 
are hi = Qu. We distinguish the diffusive observables (for which b\ ^ 0, these observables diffuse as the polymer 
center-of-mass at large times) from the non-diffusive observables (for which b\ — 0). Finally, in all the supplementary 
information, we choose the units such that D = 1, £ = 1 and k = 1. The unit of energy is fcgT = 1, the unit of time 
is C/fc, and the unit of length is the effective bond length l = ^JkjfTjk = 1. In these units, the effective Kuhn length 

is /Kuhn = a/3- 



2. Non Markovian theory in dimension d — 1 

a. Self- consistent equations for the non Markovian theory in dimension d = 1 



We now describe how to derive the equations (5,6,7) of the main text that define the non Markovian theory in 
dimension d = 1. The observable is then noted X, and we calculate the mean first passage time for X to reach the 
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value X = 0. The starting point is Eq. (3) (in the main text), which can be written in terms of modes: 



rP sta t(|a)) 



dt[P{\a),t\ir,Q) - P(|a),f|ini,0)] 



(A6) 



This equation is valid only for the modes |o) such that (b\a) — 0. We note the mathematical trick: P stSi t(\a))S({b\a) 
X) = P s tat(^)Pstat(|a)|(&|a) = X). Therefore, multiplying Eq. ^ by S((b\a)) leads to a reinterpretation: 



(0)P s tat(|a)|0) 



dt[P{0, t|?r, 0)P(\a),t\0, t; tt, 0) - P(0, i|ini, 0)P(\a),t\0, t; ini, 0)] 



(A7) 



where P(\a),t\0, t; tt, 0) is the probability of observing the configuration \a) at t given that the value of the observable 
X is X = at t and that the distribution of modes at t = was tt. Similarly, P s tat(|a)|0) is the stationary probability 
to observe a configuration given that the value of the observable is X — 0. Now, the equation (A7| is valid for any 
value of | a) (not only for those that satisfy (b\a) = 0). Noting that the distribution P(\a),t\X,t;TT,0) is normalized 
to 1, the integration of Eq. (A7) over all the modes leads to: 



TPstat(O) 



dt[P(0,t\ir,0) - P(0,t|ini,0)] 



(A8) 



Then, multiplying Eq. (A7| by and integrating over all the modes leads to: 



dt 



P(0,t|7r,0)/^'°-P(0,t|ini,0)/z£ 



(A9) 



where fxj' is the mean value of aj at t given that X = at t and that the initial distribution at t — is the splitting 
distribution tt. Similarly, multiplying Eq. (A7) by a^aj, integrating it over all the modes and using Eq. (A8) leads to 
a second set of self consistent equations: 



dt 



P(0,t|7r,0) (7 y ' +fi i ' Hj -<?ij 



r>/r\ xl' • n\ ( ini,* . ini.O ini.O 

P(0,t|mi,0) -a 



stat , * 
ij 



= 



(A10) 



where 7^' is the covariance between etj and a.j at t given that X = at t and that the initial distribution at t = is 
the splitting distribution tt. We now derive "propagation" and "projection" formulas that will be useful to explicitly 
write all the terms appearing in Eqs. ( A8|A9|A10 ). 



b. Propagation and projection formulas 



The Fokker-Planck equation that governs the evolution of the Rouse modes (in one dimension) is: 

N 



dP(\a),t) 
dt 



E 

»=i 



d_ 
da. 



AifljP- 



d_ 
dai 



P 



(All) 



It is well known that Eq. (All I admits Gaussian solutions [40] which are characterized by the average \ii of each 
mode dj and the covariance matrix 7^ that describes the correlations between the modes dj and cij . The evolution of 
/Xi and jij satisfies the following equations |40] that are sometimes called generalized fluctuation-dissipation relations 
(see [U): 



-AiUi ; %j = -(A,; + Xj)^ij + 2Sij 



(A12) 



If the initial condition is a gaussian distribution with moments and 0^ , then the value of /ii and 7^ is the solution 
of Eq. (A12) with initial conditions /i;(0) = mi and Jij(0) = cr^j. We find: 



Sn (1 



-2\it\ 



/Xi 



-\it -\jt 



(A13) 
(A14) 



These formulas describe how the mean vector and the covariance matrix are modified with time, and we call them 

2t (with Ai = 0). 



"propagation formulas". Note that Eq. (A14| is written with the convention that (1 



-2Ait 



)/Ai 
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Let us now assume that P(|a)) is a gaussian distribution, with means fii and covariance matrix 7^. We describe how 
to obtain the distribution P(\a)\X), which represents the distribution of modes given that the value of the observable 
is X. Noting that P(X\ \a)) = S((b\a) — X), and using the Bayes formula, we obtain: 



P(\a)\X)=5((b\a)-X)P(\a))/P(X) 



(A15) 



We note that the distribution P(X) is gaussian, with means (b\a) and variance (6|7|a). The distribution P(\a)\X) is 
also a gaussian distribution, and we call fj, x and j x its vector and covariance matrix. By definition, we can identify 
the value of \x x by writing: 



d\a) a z P(\a)\X) 



h\b) 



(b\7\b) 



im-x) 



(M6) 

(A17) 



The p assage from Eq 
(A15), and |e,-) represents the 1 



(A16l to (A17) results from the explicit calculation of the integral (A16l with the use of Eq. 

th 



basis vector (all its elements are except for the i which takes the value 1). The 



elements of the covariance matrix can be identified with the same method: 



y x = I d\a) (a, 



Hf)( aj -(if) P 



\X) 



Hi 



<6| 7 |6} 



(A18) 
(A19) 



Note that 7 does not depend on the value of X, which is why we ju st no te it 7 * and not 7 , to the difference of 



fj, which depends linearly on the value of X. We call the equations (A17) and (A19) "projection formulas": they 



describe how the mean and the covariance of the are modified when one restricts the modes to be on the hyperplane 
of equation (b\a) = X. We can easily see that the average of X over the distribution P(\a)\X) is (b\fi X ) = X, and 
that the covariance of X vanishes (because 7*16} = |0}), which means that X is known with certainty to be X on this 
distribution, as expected from the definition (A15). 



c. Explicit expressions of all the terms appearing in the self-consistent equations 

We remind that mj and cr^ are the mean vector and covariance matrix of the splitting distribution 7r(ai, ...,ajsr)- 
The avera ge and covariance 7^ of the modes at at t starting from 7r at t = are deduced from the "propagation 
formulas" ( |A13|A14| ): 



Mi 



< e~ Xit ; 7 ? 



%(1 



-2Ajt 



)/X i +e- x * t e-^ t af i 



(A20) 



Noting that (6|/i 7r } and (b\^\b) are the average value of X and the covariance of X at t starting from the initial 
distribution tt, we can write explicitly the value of P(X, t\ir, 0): 



P(X,t\w,0) 



{x-(b\^)y 

e 2<b| 7 T|i,> 

[27r(b\r\b)} 1/2 



(A21) 



Then, the average and covariance 7^'* of the modes at t given that one observed X at t and that the initial 
distribution is 7r are obtained by the projection formulas ( A17|A19 ): 

WI6} 



(b\r\b) 



i(b\»l-X) ; H 



(ej\^\b)(b\^\ej 
(b\T\b) 



(A22) 



For the moment, we choose gaussian initial conditions Pj n i, which has means m™ 11 and covariance matrix a. 



mi 



The moments mf 11 and a]" 1 can be chosen so that the initial distribution is the stationary distribution restricted to 
configurations such that X = X . In the case of a non-diffusive variable, this choice of initial conditions is obtained 
by applying the projection formulas ( A17|A19 1 to the moments of the stationary distribution (given by m| tat = and 
oT = <WAi): 



X bj 
XL* 



ini stat,* 



bjbj 
X l X j L 2 



(A23) 
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where we have noted L the stationary value of the root mean square of the observable X: 

N 

L 2 = {b\ a ^\b)=Y. h V\ (A24) 

In the case of a diffusive variable, one has to take care of the mode with v anishi ng eigenvalue Ai = 0. In this case, 
the moments m™ and cr™ are found by taking the limit of small Ai in Eq. (A23), we find: 



ini ini _ Mal.:i 

m l - — ; c 

ini.X 



Sij/Xi ifi,j>2 
< = = \ if 3 > 2,i = 1 (A25) 



The quantities /i i ' and 7^ - ' are deduced from to 1 " 1 and <r™ by replacing the superscript "7r" by "ini" in Eqs. 
( |A20|A22| ). Let us now define the functions (f> and tp such that: 

(X(t))™ = (% ini ) = X cf>(t) (A26) 
(AI^^-Hilfitl^W (A27) 

The function <p describes how the average of X evolves with the time when the initial condition is stationary with 
a value X = Xq, whereas ip(t) is the variance of X (with the same initial conditions). Using the projection and 
propagation formulas, we can find the values of (f> and ip. In the case of a non-diffusive variable, we obtain: 

N , 2 -Ait 

<j)(t) = V \ 6 ' ; i/j(t) = L 2 [l - 4> 2 (t)] (Non-diffusive variable) (A28) 

i=1 1 

whereas in the case of a diffusive variable, we have: 

4>{t) = 1 ; ij){t) = 2b\t + &|(1 - e~ Ajt )/A J (Diffusive variable) (A29) 

The equations ( |A28|A29 1 mean that the average value of X remains constant with time if X is diffusive, but decreases 
to a stationary value when X is not diffusive. At long times, ip ~ L 2 for a non-diffusive variable (that relaxes to a 
stationary value), whereas ip ~ &ii for a diffusive variable (that diffuses at long times). Let us finally write the explicit 
expression for P(X, t\'mi, 0): 

_ (X-Xp^,) 2 

P(X,t\iui,0)= e {2 ^ 1/2 (A30) 

All the terms appearing in the self-consistent equations (5,6,7) of the non Markovian theory are explicitly written in 
this section in Eqs. ( A21|A22|A23|A25|A30 ) . Note that the means m* associated to the positions in the main text 
are easily calculated from the means associated to the modes ai derived in this Supplementary Information (that are 
abusively represented with the same notation mf) by using the rotation matrix Q defined in Eq. (A3 1. 



3. Self-consistent equations for the non-markovian theory in dimension d = 3 

Here, we briefly describe how the theory can be extended to the 3-dimensional case. Let us assume that the 
initial conditions arc isotropic. Then, the reaction can take place anywhere on the target, and we can introduce the 
probability 7To(|a)) of reacting with a configuration |a) given that the observable has angular coordinates ft = (0, 4>) 
when the reaction takes place. Let us introduce the usual basis of unit vectors for the spherical coordinates (u r , u$, u^), 
and the coordinates (cti, r , a^e, a-i.<f>) of en in this basis: ~ ai, r u r + ai^giig + a^u^. Due to the symmetry, we can 
separate the coordinates in the splitting distribution: 

7ro(|a)) = 7r I .(|a I .))7r e (|a e ))7T0(|a0)) (A31) 

The gaussian approximation is written for each coordinate: we assume that the radial splitting distribution n r is a 
multivariate gaussian, with mean vector and a covariance matrix cr^ , whereas the two perpendicular splitting 
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distributions irg and ir^, are multivariate gaussian, with vanishing mean vector (mj' — 0), and a covariance matrix 
crjj . The resulting self-consistent equations are very cumbersome, that is why in this paper we make the simplifying 
assumption of isotropy of the covariance matrix: we assume that a 71 '" = <t' t '- l , and we note the common value of these 
matrices a™ . We note \i\ and 7^ the radial mean vector and the covariance matrix of the distribution P(|a), i|7TQ, 0). 
Let us specifically assume that the initial distribution is the stationary distribution restricted to configurations such 
that ||R|| = Rq. Let P(|a) , i|Pi n i,n, 0) be the distribution of modes at t starting from an initial stationary distribution 
where the observable value is Ro = RoU r (8, <j>) and has the angular spherical coordinates = (8, (j>). As in Eq. ( A31 1, 
we can separate this distribution into P r PgP<j,, where each function is a gaussian. We now rewrite the equation (4) in 
the main text in the space of modes and with an explicit average over the angles: 

rP stat (|a)) = y dtJdSl[P(\a),t\ir n ,0)-P(\a),t\IU,a,0)] (A32) 
where dQ = sin 8 d8d(f>/ (Air). We multiply both members by <5((c|a) — Rf) and reinterpret this equation: 

rP stat (Rf)P sta t(|a> |R f ) = J dtj dO[P(R f , t|7r n , 0)P(|a), <|R f , t; 7r , 0) - P(R f , i|ini, 0)P(|a), t\R f ,t; P ini , n , 0)] 

(A33) 



This equation is the equivalent in 3 dimensions of Eq. ( A7 ) . Integrating it over all the modes gives an estimation of 
the mean first reaction time: 



rP stat (R f ) = J dtj dn[P(Rf,t|7r n ,0) - P(R f , t|P mi ,n, 0)] 



(A34) 



For symmetry reasons, we can assume that Rf = ZfU z without loss of generality. The distribution P(Rf, t|7rn, 0) can 
be calculated by noting that for a given value of 8, the radial component of Rf is Z{ cos 8 whereas its component in 
the 6 direction is — Z{ sind: 



P(Z f u,, i|7rn, 0) = P(cos 0Z f , t\ir r , 0)P(- sin 6Z { , t\n 9 , 0)P(0, t\n^, 0) 



1 (Z f co39-(6| M ' r >)' ! , (-Z 



In the same way, we have: 



P(ZfU„t|P mi ,n,0) 



1 

[2ir(b\~/"\b)] 3 / 2 

e - i [(Z f co S 9-i?. o 0) 2 + (Z fS in9) 2 



[2^1 2 



Obviously, the equation (A34) can be simplified by taking Zf — 0, in which case we obtain: 

poo 

rP sta t(0)= / tft[P(0,t|7r n ,0)-P(0,t|fl n i 1 n,0)] 



with: 



P(0,t|7r n ,0) = 
P(0,t\PU, n ,0) 



1 



[2ir(b\j*\b)} 3 / 2 

1 Cgoj) 



(A35) 
(A36) 

(A37) 

(A38) 

(A39) 
(A40) 



[2tt^] 3 / 2 

The self-consistent equations that define the values of ni* are obtained by multiplying Eq. (A33) by a^ 2 and by 
integrating over all the modes. To do so, we first need to evaluate the integral: 

J d\&)a iz P(\a),t\Z{U z ,t;iT n ,0) = J d\a)(a lr cos8 - a i8 sin0)P(|a), t\Z { u z , t; Tr n , 0) 

W\b) 



cosdrfl^ + Z[- 



■c\r\b) 



(A41) 



In the passage from the first to the second line, we used the fact that Oi Z = a ir cos8 — etjgsin^. The passage to the 
third line uses two times the projection formula ( A17). In principle, we can choose any value for Zf € [0; a]. In order 
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to obtain simple formulas, we chose to develop the equations for the first moment at lowest order in Zf. Reporting 
Eq. (A41) into the equation (A33l (multiplied by a,- z ) and integrating over all the modes leads at lowest order in Zf 
to: 



rP S (0)Ki 



dt 



1 o (ejlYW 

(p) + M l?(Mk ' 0) 



(A42) 



where 



Ki = 



$ii/bi (diffusive variable) 
bi/(XiL 2 ) (non-diffusive variable) 



(A43) 



We do not give any details on the derivation of the equations that define the second moments at the limit Zf = 0, we 
multiply Eq. (A33) by ai Z aj z and integrate over the modes, we obtain: 



rP stat (0)<7: 



stat,* 



(it 



' 7T,0 7T,0 \ 

^-^ + 7 £* P(0,f|7r„,0) 



' ini,0 ini,0 
IH mi,* 

3 



iij 



P(0,t\P ini , n ,0) 



(A44) 



The equations (A38), (A42) and (A44) form a set of non-linear equations that enable to compute both the mean 



first reaction time and the moments of the splitting probability distribution under the hypothesis of isotropy of the 
covariance matrix. They are the generalization to a 3-dimensional space of the equations (5,6,7) of the main text. 



4. Mean first passage time averaged over initial conditions (non-diffusive observable) 



In the case of a non-diffusive variable, we can ask for the time it takes for the reaction when the initial distribution 
is the stationary distribution (restricted to configurations that lie outside the reactive region, ||R|| > a). Then, the 
initial distribution is a superposition of the distributions restricted to a given value of Rq: 



POO 

P mi (\a)) = Prtat(|o)|||R|| > o) = / dRoK%P stat (Ro)P s t at (\a)\\\IL\\ = R ) 

J a 



with: 



e -R 2 /(2L 2 ) 

Pstat(Po)= z{a L2) ; Z(a,h) 



dR R 2 e- R WV 



(A45) 



(A46) 



Eq. ( A45| states that the initial distribution is a superposition of the initial distributions studied in the section A3 



We can therefore follow step by step the calculations that led to the equations ( A38 1 , ( A42 ) and ( A44 ) and average 
over Rq at the end of the calculation. For simplicity, we only give the self-consistent equations for the first moments 



of the splitting when the covariance matrix crfj is approximated by u i 



This approximation turns out to be an 



excellent approximation. Under this approximation, the equation that defines the mean vector of the splitting is: 

bi 



dt< 



1 «,o(bm , (eib {Xo ' ini} \b) 
3 Ml V V> 



3i/> 



he 



-Ait 



XiL 2 

e .| 7 {X ,ini}|^\ Q( a ^ 



P(0,t|7r n ,0) 

e .| 7 {X ,mi}| 6 ) 



Z(a,L 2 ) 



Z{a>,j>) 
XiL 2 J Z(a,L 2 ) 



1 



(2tt^) 3 / 2 



= 
(A47) 



where the superscript u {Rq, ini}" simply indicates that one refers to the only initial conformations such that ||R|| = Rq, 
the function Z is defined in Eq. ( |A46 1 and G is defined by: 

(A48) 



G(a, h)= dR Q R 4 e~ B l/W = a h(a 2 + Zh^/^ + iJ*-h b ' 2 [l - Erf (a/Vth) 
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The mean first passage time can be evaluated by using the formula: 



rPat,t(Rf) 
or, taking Rf = 0: 



dt / dD, 



P(Rf,t|7r n ,0) 



dR R 2 P stat {R )P(-R l ,t\{R ,im},n,t = 0) 



rP 8t at(0)= / 
Jo 



dt 



P(O,t|7r n ,0) 



Z(a, tp) 



(A49) 



(A50) 



Z(a,L 2 )(2TT^) 3 /\ 

We note that, in the limit a -C L, we have Z(a,ip)/Z(a, L 2 ) ~ tp 3 / 2 /L 3 and we obtain the simpler expression from 
Eq. ((A49b: 



/•OO 

rPstat(Rf) = / d* [P(Rf, tkn, 0) - P sta t(R f )] 
Jo 



(A51) 



This expression enables us to make the link between the Markovian theory and the Wilemski-Fixman theory with 



a delta-sink function. The Markovian approximation is written as: <7^ = a 



stat,* 



and 



stat,a 



Reporting these 



approximations into Eq. (A51) and taking Zf = a leads to the expression of the Markovian approximation of r in 
the case of a small target: 



dt 



L 2 



i^e sinh 



a 2 ^ 



L 2 {l~<j) 2 ) 



- 1 



(A52) 



This expression is equivalent to the Wilemski-Fixman theory with the delta-sink approximation introduced in Ref. 
[29]. This shows that the Markovian approximation leads to the same results as the Wilemski-Fixman theory. 



Appendix B: Scaling relations 

In this section, we derive scaling relations in 3 dimensions. We always derive scaling relations by assuming that 
the second moment of the splitting distribution can be approximated by its stationary value (cr 7r ~ cr 8 * 8 *'*). We first 
identify the relevant regimes by looking at the time scales. 



1. Properties of the walk at different time scales 

At short times, in both diffusive and non-diffusive cases, we have: 

i>(t) ~ 2D t ; D = (b\b) (t -> 0) (Bl) 

Hence, the walk is always diffusive at short times with an effective diffusion coefficient Dq. In the case that the 
observable is the first monomer position, Do = D, whereas Dq = 2D in the case of the end-to-end vector. At long 
times, for a non-diffusive observable, we have: 

ij-L 2 ; <P^b 2 2 e- X2t /(L 2 \ 2 ) (B2) 

This clearly shows that the observable relaxes to a stationary distribution. In the case of a diffusive variable, we have: 

ijj(t) = 2b\t (B3) 

In this case, the motion is diffusive at long times. In the case that the observable is the first monomer position, the 
diffusion coefficient at large times is simply the diffusion coefficient of the center of mass b\ = £>cm = D/N . Irregular 
behavior of i\> and 4> appears at the limit of infinite N. In this case, one sums over an infinite number of modes, and 



the eigenvalues given by Eq. ( A2 1 become in the large N limit 



A, ~ (g - l)7r R (B4) 

where tr = N 2 /tt 2 is the Rouse time, the slowest relaxation time of the non-diffusive modes. In the case of the 
position of the first monomer, b\ — 1/s/N and &i>2 — \/2/N. Let us call -0oo and <^oo the functions ip and <p obtained 
in the limit of an infinite number of modes. The short time limit of V'oo can be evaluated by transforming the (infinite) 
sums ( |A28|A29[ ) into integrals (see [3T] for details): 



This behavior defines an effective walk dimension: tp 




(B5) 

ri 
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2. The limit of small target size a (at fixed N) 

a. Case of a non- diffusive variable 
Let us write Eq. (A47l in the case that a n ~ cr stat '* and in the case that a -C L. 

r&i(l - e"^ 



r 

3^ 



tpXi 



A, 



1 - e 



-Ait, 



1 

L2 



e 2 * 
(2ttV') 3 / 2 



= 



(B6) 



where we have set T = (b\fj, w ) and we have used the fact that: -yl stat >*} = j{X ,uu} anc i that: (ei|7^ stat '*^ |6) = 
bi/\i(l — e~ Xit (f>). We assume the scaling = a fhi ; T(t) = a T(t) for a — >• 0. We check this scaling by calculating 
all the terms of Eq. (B6| in the limit a —¥ by taking their short time expression. For example, we have: 



00 Tmf r e~ Xit 
dt l ' r 



Zip (27n/>) 3/2 



dt 



a 2 fhi e 



-a 2 /[2(2D Q t)] 



3(2D Q t) [27r(2L> t)] 3 / 2 12irD a 



(B7) 



All the other terms of Eq. (B6) are also of order 0(a 1 ), which proves that the development for small a is consistent. 
Equating all the coefficients of a -1 gives the following estimation for the moments: 



a mi 



A; 



A^ 
A) 



l 

Z 2 " 



(B8) 



Note that this formula is consistent with the condition (6|m lr ) = a (because (b\b) = Do)- The conclusion is that the 
moments of the splitting probability vanish for small target size, and therefore they do not play any role in the global 
mean reaction time: the Markovian theory gives the same result as the non Markovian theory in this limit. This time 
can be evaluated by replacing mf by in Eq. (A50) and then evaluating for small a the resulting integral: 



rP«tat(0) 



dt 



e -a 2 /[2(2D t)] 

[2n(2D t)}3/z 



1 



4irDna 



(B9) 



This behavior of r is the same in the Markovian theory and the non Markovian theory and is fully compatible with 
the results of the other classical Markovian theory (the harmonic spring model [3l[25]). The monomers behave as if 
they were disconnected in an effective volume V^f = 1/-P s tat(0). 



b. Case of a diffusive variable 

We now derive the behavior of t and in the case of a diffusive variable in the limit a — > at fixed initial distance 
Rq. We write Eq. (A42| for i > 2 for a large initial distance Ro — > oo: 



dt 



r 



r&i(i 



-AiiV 



IpXi 



e 2 * 

(2^)3/2 



(BIO) 



where we have set T = (blfi*). We apply again the method that we described in the case of a non-diffusive observable. 
We assume the scaling mf = a fhi, T(t) — a T(t) for a 0. This scaling is consistent and leads to the following 
evaluation of the moments: 



-ah/ D if i > 2 

a{2/b 1 -b 1 /D ) if i = 1 



(Bll) 



where the second equality is deduced from the condition (6|m 7r ) = a. Note that, in the case that the observable is the 
first monomer, this equation implies that the average position of the center-of-mass at the reaction is: 



(#cm)7r = b x mJ[ = a 1 + 



,2 

i=2 °i 



bi + Ef= 2 bi 



> a 



(B12) 
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That is to say, the position of the center of mass at the reaction is at the exterior of the target, at a distance of order 
a, as expected from the intuition. In this case also, the average reaction time is evaluated to be: 

f°° e ~ o 2 / [2(2 D t)] y 

T ^ V Jo A [2«{2D t)]W = 4^a (B13) 

This last equation shows that, in this regime, the first monomer behaves as if it were not connected to the rest of the 
polymer chain. 



3. Scaling in the thermodynamic limit TV — > oo 



Reaction between the first monomer and a target. 



Here, we determine the scaling relations for the mean first passage time at the thermodynamic limit TV — > oo in 
the case that the observable is R = rj. In the limit of large TV, the number of modes that must be taken account is 
infinite, the Rouse eigenvalues are approximated by A g ~ (q — l) 2 ir 2 /N 2 , and the coefficients b q are approximated by 
b q > 2 ~ \j2jN. We introduce the rescaled variables rh q = m q /N, a = a/y/N, f = t/N 2 , b q = b q VN, V = V/N 3 / 2 . In 



\J 2 1 TV. HV I L I I ■ I ' " L L V. ' I 1 ! V .M (ill 11 w I J. I (II ' n ^ in,, 

terms of these rescaled variables, one obtains the equations valid for q > 1 



dt 



r 

3-00 



iq+ie 



f%/2(l - 



ooq 2 n 2 



v/2(l - e-g 2 " 2f ) 
^oog 2 7r 2 



the first moment rh\ is determined such that fn\ + V2^2 q>2 ~ S, which leads to 



e 2 </-o 



(27T0OO) 3 / 1 



= 



(B14) 



°° °° i _ P -q 2 n 2 t 

f = 5 + ^m- +1 (e- w qt -l); oo = 2i + 4^ 



9=1 



9=1 



and the evaluation of the mean first passage time is 



q 2 Tr 2 



(B15) 



(B16) 



This defines a scaling relation t/V — yNf(d) = VNf(a/VN), where / is a dimensionless fun ction that can be 
evaluated numerically. The first numerical approach consists in trying to solve the exact equations ( A42 1 for TV large. 



This method gives the exact value of the mean first reaction time for any finite TV, but is not suitable to determine 
the mean first reaction time in the limit of a small target. Indeed, it can be seen that all the curves obtained with 
this method converge when TV oo to a single curve /(a) (Fig pi), but that the convergence is very slow for small 
target sizes. This method is therefore not suitable to estimate f[a) for 5 < 0.2. The reason of this difficulty comes 
from the finite size effect detailed in section B 2 in the limit of small size a -C 1/TV, the mean first passage time scales 
as 1/5. To overcome this difficulty, we directly approach the solution of the rescaled equations (B14| by introducing 
a cutoff TV C beyond which we approximate the moments by (m q = if q > N c ). The solutions obtained for finite TV C 
are expected to converge for TV C — ¥ oo, leading to the determination of f{a) without having the problem of the finite 
size effect, as can be seen on Fig [6} This method enables us to determine the whole function /(a). 

The scaling in the Markovian approach is much more straightforward. Inserting Y ~ a into Eq. (B16) yields: 



e 2 ^oo 



Vfi 



Markovian 



(a) 



(B17) 



Therefore, the Markovian estimate for /(O) is /(O) = J dt(2mp 00 )~ 3 ^ 2 ~ 0.112, which gives the estimate a e ff/\/TV = 
l/[47r/(0)] = 0.71. The difference with the non Markovian result is approximately 21%, and reaches 100% for 
intermediate target sizes. The effective target size is related to the function / by: 



a c ff 



TV 
47r/(a) 



(B18) 
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FIG. 6: Determination of the scaling relation for the intermolecular reaction time averaged over initial distances. 

Circles: solution of Eq. ( |A42[ ) valid for finite iV for various values of N (represented with a color code). The dashed lines 
represent the prediction of the Markovian theory for the same values of parameters. The black stars represent the solution 
of Eq. ( B19| l when truncated at a large enough N c and give an estimation of the function / [Eq. ( |B16| )] , while /Markovian is 
represented by the black dashed line. 



b. Cyclization reaction 

In the case of the cyclization reaction, we apply the same method. The rescaled equations read, for q odd: 



dt 



-q 7T t. 



q 2 ir 2 



4>o 



e 2 </>o 



(27rVoo) 3 / 2 



(B19) 



with: 



q— l,odd q— l.odd 

and the evaluation of the mean first passage time is : 



q 2 it 2 



(B20) 



t —- (2tt) 3 / 2 /"dt I 6 " 2 "°° 1 1 - ^ 



(B21) 



We apply the same method as in the case of the intermolecular reaction to determine the dimensionlcss function c. 
One additional difficulty in this case is that the truncated equations do not allow a solution such the solution such that 



? =1 



a because the redundancy of the equations is lost when the equations are truncated. We have therefore 



E 

to release one equation and to impose the condition E 9 =i m q ~ ^- We tried to leave the restriction on the first mode 
as well as the last one, leading to almost no difference in the limit of large N. The results are presented on Fig. [7J 
We find the scaling r ~ 1.732r^ in the limit of small target size. The Markovian estimation gives: 



^Markovian — — (271") 



3/2 



dt 



(27n/>oo) 3 / 2 (2tt)3/2 



^Markovian (a) 

IT 2 



(B22) 



From Eq. (B22), we estimate the Markovian scaling relation for very small target sizes to be r ~ 1.977tr for small 



target sizes, the difference with the non Markovian scaling is about 14%. For realistic target sizes (for example for 
a = 0.2), c ~ 0.5cMarkovian: the difference between the Markovian and non Markovian scaling is of the order of 100%. 
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FIG. 7: Determination of the scaling relation for time the cyclization time averaged over stationary initial 
conditions. Continuous lines: global mean cyclization time predicted by the non Markovian theory, obtained by solving Eq. 
(A42| for various values of N (represented with a color code). The dashed lines represent the prediction of the Markovian 

The black stars represent the solution of Eq. (B19l when truncated at a large enough N c , 
The black dashed line represents CMarkovian [Eq. (B22l]. 



theory for the same values of N. 
and represent the function c appearing in Eq. 



IB21I 



Appendix C: Details of simulations 
1. Simulation of cyclization reaction 



We performed simulations of cyclization events by using the algorithm presented in Ref . |29j . The initial condition 
is an equilibrium configuration of the polymer, in which the coordinates of each mode a,i tX i a >i,y> a >i,z can be chosen 
from a normal distribution with variance 1 /\ and zero mean. The initial positions of the monomers are deduced from 
the initial values of the modes by using the transformation rule of Eq. ( A4 1 . If the initial configuration is inside the 
reactive zone (|r ce | < a), this configuration is rejected and another configuration is chosen. After the determination 
of the initial configuration, the monomer positions evolves at each time step with the following algorithm: 



Xi{t k +l) = Xi(tk) - 



N 



(At) k + y/2 (Ai) fc 



(CI) 



where )X is a number generated with the normal distribution with variance unity. The evolution of the other 
coordinates follows the same equation, and the time step is variable [29 : 



{At) k — £ feH 



(At) low + (At) high sin (|R ee (i fc )| 2 - a 2 ) tt/6) 
(At)i ow + (Ai) high 



if |R ee (t fc )| 2 <a 2 + 3 
otherwise 



(C2) 



With this choice, the time steps become smaller and smaller when the reactive zone is approached. The slight 
differences of numerical constants with Ref. [29] is due to a slightly different choice of units. The simulation runs 
until the condition |R ee | < a is satisfied, in which case the value of r = ^2u(At)k i s an estimation of the first cyclization 
time for this trajectory. We have chosen (Ai)hi g h = 10 -4 , and (At)i ow = 10 -7 . These values are approximately 3 
times smaller than the corresponding values in Ref. |29j . 



2. Simulation of intermolecular reactions 



We now describe the simulations that lead to the estimation of the mean first passage time of the first monomer 
to a target located at x = 0, when there is a reflecting wall at x = L (we first expose the method in d = 1). In the 
simulations, only the first monomer is affected by the presence of the target and of the reflecting walls. The initial 
value of a mode (i > 2) is taken from a normal distribution of variance 1/A,. The initial value of a\ is chosen such 
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that the position of the first monomer is x\. The initial position of the monomers is then obtained by applying Eq. 
(A4). The, the position of the monomers Xi,i > 2 evolves at each time step with: 



Xi(t 



k+l) 



Xi{tk) 



N 

E 



MijXj(t k ) 



(At) + y/2 (At) u lM 



(C3) 



in which u^ k is a random number taken from a gaussian distribution with variance 1. If Xx(t k ) is far from the 
boundaries, it also evolves according to this equation: 



JY 



xi(t fc+ i) = xx(t k ) + FxAt + V2At ui, k ; Fx = -^AfySjftfc) 



(C4) 



If the position x\ is close from the reflecting wall, then the last equation has to be modified in order to take into 
account the fact that the particle has a decreased probability to approach the wall (and zero probability to cross it). 
Therefore, if \L — x±\ < d, we follow the procedure introduced in Ref. [55] : 



Xl (t k+1 ) = X X (t k ) - /{' cfl ( ^yp-\ ^At + Ui, k / 2 rcfl (L " 



At + F x At 



(C5) 



In this equation, ux.k is a random number that takes the values ±1 with equal probability, and the positive functions 
/[ efl and /j 03 are the functions fx and / 2 of the equation (18) in the reference [3B]. Similarly, if Xl is close from 
the absorbing boundary (i.e. xx < d), one calculates P a b s = 1 — ert(xx/ (2V~Ai)) the probability of being absorbed 
between t and t + At. One then generates a random number between and 1 to decide whether or not the target is 
reached during the time step, in which case the simulation stops. If the absorbing wall is not reached, then xx evolves 
according to: 



xx(t k + 1) = xx(t k ) + f{ 



abs 



Xx 



At + Ul , k /£ 



abs 



Xx 



At + Fx At 



(C6) 



where the random number u lyk takes again the values ±1 with equal probability. The positive functions /f bs and /.f bs 
are the functions fx and / 2 of the equation (16) in the reference [36 . At the end of the simulation, the positions of the 
monomers are recorded, thereby giving an access to the splitting probability. The parameters used in the simulations 
are: d = L/10, At = 5 x 1CT 4 , L = 158 (the volume is V = 2L = 316), N = 10. 

In 3 dimensions, we apply the same algorithm except that the equations ( C5|C6 ) are applied to the radial component 
of the position of the first monomer. The target is a sphere of radius a = 1.7 at the center of the volume (sphere of 
radius R — 28.4). The time step for these simulations is also At = 5 x 10~ 4 , and the distance d is 2.85. In all the 
simulations, the random numbers are generated with the "ran2" algorithm described in the Numerical Recipes in C 
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